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ABSTRACT 

We present magnetic flux measurements in seven rapidly rotating M dwarfs. Our sample stars have 
X-ray and H-alpha emission indicative of saturated emission, i.e., emission at a high level independent 
of rotation rate. Our measurements are made using near-infrared FeH molecular spectra observed 
with HIRES at Keck. Because of their large convective overturn times, the rotation velocity of M 
stars with small Rossby numbers is relatively slow and does not hamper the measurement of Zeeman 
splitting. The Rossby numbers of our sample stars are as small as 0.01. All our sample stars exhibit 
magnetic flux of kilo-Gauss strength. We find that the magnetic flux saturates in the same regime 
as saturation of coronal and chromospheric emission, at a critical Rossby number of around 0.1. The 
filling factors of both field and emission are near unity by then. We conclude that the strength of 
surface magnetic fields remains independent of rotation rate below that; making the Rossby number 
yet smaller by a factor of ten has little effect. These saturated M-star dynamos generate an integrated 
magnetic flux of roughly 3 kC, with a scatter of about 1 kG. The relation between emission and flux 
also has substantial scatter. 

Subject headings: stars: activity — stars: late-type — stars: magnetic fields — stars: rotation 



1. INTRODUCTION 

Stellar magnetic activity is closely related to rota- 
tion in solar-type stars. In slow rotators, activity scales 
with the rotation rate until it becomes saturated at 
a certain velocity, which means tha t it does not grow 
further regardless of rotation rate (|Noves et al.l 119841 : 
iPizzolato et"alll2003[) . This effect is seen in a broad 
variety of activity indicators (|Vilhul 11984 Il987f) . The 
threshold rate at which saturation occurs depends on 
the spectral type of the star, with the convective over- 
turning time perhaps determining this threshold. Sat- 
uration sets in where the value of the Rossby number 
Ro = P/tcouv ^ 0.1, i.e. where the timescale of rota- 
tion is significantly shorter than the timescale of typical 
convective eddies. 

Indicators of stellar activity are usually coronal or chro- 
mospheric emission observed at X-ray, UV, optical, in- 
frared, or radio wavelengths. We know from the Sun that 
this emission is induced by magnetic fields heating the 
upper layers of the solar atmosphere, and by analogy we 
conclude that stellar activity is connected to magnetic 
fields on the surface of other stars. At high rotation 
rates (or small Rossby numbers), all activity indicators 
saturate, i.e. they do not grow over a cer tain level re- 
gardl ess of higher rotation rates (see, e.g., iJames et al.l 
l2000t) . Two possible explanations exist for the satura- 
tion: 1) The stellar dynamo process saturates and stars 
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cannot produce magnetic fields stronger than the satu- 
ration value; or 2) The magnetic fields continue to grow 
at more rapid rotation, but the fraction of the surface 
filled with fields - or the area covered by spots - reaches 
unity so that no more emitting plasma can be placed on 
the star. The only way to decide which way the stars go 
is to directly measure the magnetic field. Unfortunately, 
this is very difficult and the picture - particularly in stars 
w ith sa t urate d activity - is not yet clear. 

ISaad ()1996f ) has collected measurements of magnetic 
flux Bf, of the filling factor /, and of rotation periods 
to investigate the behavior of magnetic fields on stars. 
It is important to realize that the magnetic flux, or the 
magnetic field average over the whole surface, is not the 
same as the local field. It is the (unsigned) mean av- 
erage of the magnetic field strength over the whole sur- 
face. Furthermore, the strongest magnetic fields in cool 
spots may not be fully captured because their contribu- 
tion to the total flux is dimished due to their low temper- 
ature. As on the Sun, magnetic flux is probably concen- 
trate d in relativ ely small regio ns of strong fields (see, e.g. 
John s-Krull fc Valcnti 2000). ISaad ([1996) shows that in 
stars rotating slower than the saturation threshold the 
magnetic flux Bf as well as the filling factor / show the 
same trend as all other activity indicators: They grow 
with larger rotation rate. At high rotation rates, ISaail 
(1996) claims that saturation occurs in the filling fac- 
tor / but not in the magnetic flux Bf . From this result 
one would conclude that a star rotating at the saturation 
threshold is completely covered with magnetism (/ = 1), 
and that the saturation phenomenon is due to the satu- 
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ration of the emission process while the sta r's m a gnetic 
flux can grow further with higher rotation. ISaad (|2001l ) 
has reinvestigated this issue with a few more data points, 
noting that there is some indication for a saturation at 
Bf ~ 3kG at small Rossby numbers. Much higher mag- 
netic flux could be in contradiction to the idea that mag- 
netic fields in stellar atmospheres cannot grow stronger 
than the equipartition field, i.e. the field strength at 
which th e magne ti c pre ssure equals the gas pressure. 
However, ISolankil ()1994) shows that the eqiupartition 
field may not necessarily be a hard upper limit for the 
field strength at / = 1, so that more rapidly rotating 
stars could in principle have much stronger fields. 

The only way to decide whether Bf does saturate or 
not is to provide direct measurements of magnetic flux in 
the regime of saturated activity. Generally, the measure- 
ment of magnetic fields relies on the splitting of spec tral 
lines through the Zeeman effect fe.g. lRobinson|[l980l ). In 
rapid rotators, the subtle effect of magnetic broadening 
is buried under the rotational line broadening so that it 
is particularly difficult to directly measure the magnetic 
flux in the regime of saturated magnetic activity. How- 
ever, the crucial datum for rotational line broadening is 
the projected rotation velocity v s'mi and not the Rossby 
number or the rotation period (which are the numbers 
that appear to set the saturation threshold of magnetic 
activity) . Mainly because of the smaller radius of cooler 
stars, the surface rotation velocity at which saturation 
sets in depends on spectral type. In early G-type stars, 
the saturation velocity in the transition region is on the 
ord er of 30kms^^ (and only 15kms~^ in the corona, 
Avrcs 1999). In M dwarfs, it is less than Skms"""^ (see e.g. 
Reinera,2007i) . Thus, in sun- like stars, the high surface 
rotation velocity required for activity saturation hampers 
the measurement of magnetic flux, but this regime can 
easily be probed in M stars. 

One potential problem with the use of M dwarfs for 
the investigation of dynamo related phenomena is that 
the interior structure changes around spectral type M3.5 
- stars cooler than that are completely convective. Nev- 
ertheless, no change in activity is observed at the thresh- 
old to complete convection. A rotation-activity connec- 
tion is observed in M stars down to sp ectral types M 8.5 
(jMohantv fc Basrill2003HReinersll2007l ). 

In this paper, we present direct measurements of mag- 
netic flux in several M stars. Some of them are very 
rapid rotators and clearly belong to the regime of satu- 
rated activity. We aim to clarify whether the magnetic 
flux Bf saturates as Ha and X-ray emission do, or if Bf 
continues to grow beyond Bf ^ 3kG. 

2. SAMPLE AND OBSERVATIONS 

For our sample we chose a number of mid-M stars 
with known X-ray emission and presumably high rota- 
tion r ates. The valu e s f sin^ and logLx/^boi are taken 
from iDelfosse et al.l except for GJ3379. This 

value is calculated from the X-ray luminosity taken from 
ISchmitt fc Liefkd (|2004f ). With one exception, projected 
rotation velocities v sini were available for all targets. 
We chose only stars in which v sini was reported to be 
above 5kms~^ and which show saturated normalized X- 
ray emission. For GJ 3379 we are not aware of any former 
V sini measurement, but the high value of normalized X- 
ray emission is indicative of saturation and we added the 



TABLE 1 
Log of observations. 



Name 


UTC Date 


Exp. Time [s] 




GJ 3379 


2007-09-30 


200 


-2.86 


GJ 2069 B 


2008-01-24 


600 


-2.77 


Gl 493.1 


2007-04-25 


600 


-3.31 


LHS 3376 


2007-04-25 


1800 


-3.63 


GJ 1154 A 


2007-04-25 


600 


-3.28 


GJ 1156 


2007-04-25 


600 


-3.39 


Gl 412 B 


2007-04-25 


1200 


-3.28 



TABLE 2 

Results of ou r analysis. Literature values are from 
IDelfosse et al.I II1998I) . 



Name 


SpT 


V sini 


Bf 




V siniLit 






[kms^-'-] 


[G] 




[kms^^] 


GJ 3379 


M3.5 


< 3 


2300 


-3.35 




GJ 2069 B 


M4.0 


6 


2700 


-3.28 


9.2 


Gl 493.1 


M4.5 


18 


2100 


-3.46 


16.8 


LHS 3376 


M4.5 


19 


2000 


-3.73 


14.6 


GJ 1154 A 


M5.0 


6 


2100 


-3.55 


5.2 


GJ 1156 


M5.0 


17 


2100 


-3.53 


6.5 


Gl 412 B 


M6.0 


5 


>3900 


-3.72 


7.7 



star to our sample. 

Data were taken at the W.M. Keck observatory with 
the HIRES spectrograph. Our setup covers the wave- 
length range from below Ha (6560 A) up to the molec- 
ular absorption band of FeH around 1 /J,m. We used a 
slit width of 1.15 "achieving a resolving power of about 
i? = 31 000. Our log of observations is given in Table[TJ 
Data were cosmic-ray corrected, flatfielded, background 
subtracted, and wavelength calibrated using a ThAr 
spectrum. Data reduction was carried out using routines 
from the echelle package within the ESO/MIDAS dis- 
tribution. Fringing is not an issue in spectra taken with 
the new HIRES CCD, even in very red spectral regions 
around 1 /im. 

3. ANALYSIS 

The analysis of our spectra follows the strategy laid 
out in lReiners fc Basril ([200611 and Reiners & Basri (2007, 
RB07 in the following). We measure the equivalent width 
of the Ha emission and convert this number to normal- 
ized Ha luminosity using M-star atmospheres calculated 
with the PHOENIX code (|Allard et all 1200 1[ ). To mea- 
sure the projected rotation velocity v sini and the mag- 
netic flux Bf of our sample stars, we utilize the absorp- 
tion band of molecular FeH close to 1 /xm. We com- 
pare our data to spectra of the slowly rotating M-stars 
GJ1002 (M5.5) and 01873^ (M3.5). In order to match 
the absorption strength of the target spectra, the inten- 
sity of the FeH absorption lines in the two comparison 
spec tra is modifi ed according to an optical-depth scaling 
(see iReiners fc Basri 200^. In a first step, we compare 
the artificially broadenend spectrum of GJ 1002 to the 
target spectra in the wavelength region at 9930-9960 A 
to determine the value of v sini by x^-minimization. 

G1873 is rotating a t v sini ~3km~^ (R B07), not at a higher 
velocity as reported in IDelfosse et all l|1998l ^. The small (but de- 
tectable) rotation of Gl 873 does not affect our measurements. 
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For the determination of the magnetic flux Bf, we 
concentrate on smaUer wavelength regions that contain 
absorption hnes particularly useful for this purpose, i.e. 
regions that contain some magnetically sensitive as well 
as magnetically insensitive lines. The magnetic flux of 
G18 73 was measured to be 3.9 k G (using an atomic Fel 
line: lJohns-Krun fc Valeirtill2000f) . For our measurement, 
we are using a spectrum that contains both the Fel line 
and the FeH absorption band. The Fel line in this spec- 
trum is consis tent with the same magnetic flux value as 
found bv Johns-Krull fc Valenti (2000) so that we can use 
the FeH pattern for the calibration of magnetic flux mea- 
surements in other stars. This method does not require 
theoretical models of the magnetic Zeeman splitting of 
FeH lines, which are not available yet. 

We determine the magnetic flux of our target stars by 
comparison of the spectral regions at 9895.5-9905.5 A, 
9937.5-9941.0A, 9946.0-9956.0A, and 9971. 5-9981. OA 
(for more details see RB07). In Figs.[TH4] we show the 
data and the quality of our fit in the top panels. Note 
that in rapid rotators the difference between magnetic 
and non-magnetic stars is not necessarily clearest at the 
exact location of magnetically sensitive lines. The blend- 
ing of lines through rotation pronounces differences at 
wavelengths where the equivalent widths of lines differ 
the most between active and inactive stars. 

For example, the two FeH lines at 9949.1 A and 
995 1.7 A are magne t ically not very sensitive (see 
iReiners fc Basril I2006L 120081 ). but it appears that the 
9949.1 A line effectively gains a little in equivalent width. 
Thus, at rotation rates as high as shown in Figs. [3] and [H 
the region around 9949 A becomes the one of largest dif- 
ference. The reason for this are the magnetically sensitive 
features next to insensitive lines together with the effec- 
tive gain in equivalent width, although such differences 
are not necessarily overt for the observer once Dopplcr 
broadening is introduced. On the other hand, it can also 
happen that at wavelengths where the unrotated spectra 
are quite different, rotational broadening averages in ad- 
jacent flux in such a way that the magnetic differences 
end up erasing each other there (e.g. at 9948 A). Thus, it 
is necessary to carry out a spectral fitting procedure af- 
ter proper preparation instead of relying on appearances 
in original template spectra at fixed wavelenghts. 

4. RESULTS 

In the bottom panels of Figs.[TH4l we show the x^- 
landscapes for all our targets as a function of v sin i and 
Bf. Color-coding displays the quality of the fit. In each 
X^-landscape, the white contour marks the 3a region, i.e. 

< Xmin + 9 within this region. In all cases, the mean 
deviation per degree of freedom is on the order of 1 (xi^ ~ 
1) for the estimated signal-to-noise ratio. Uncertainties 
in V sinz and Bf are typically around lkms~^ and a 
few hundred Gauss, respectively. We emphasize that in 
particular in the case of Bf systematic errors are a more 
severe source of uncertainty so that the total uncertainty 
in Bf is more realistically in the 500-1000 G range. 

The results of our analysis are given in Tabled pro- 
jected rotation velocity v sini, magnetic fiux or the 
mean magnetic field Bf, and normalized Ha activity 
log L-Ha/Lho\ are given in columns 3, 4, and 5, respec- 
tively. All Ha measurements confirm that our sam- 
ple targets are active stars close to the activity satu- 



ration level. Three stars show very rapid rotation on 
the order of 20kms~^, three stars are rotating at a ve- 
locity around 6kms~^. GJ 3376 shows rotation below 
our detection limit of v sini « 3kms^^. For compari- 
son, we includ e in c olumn 6 measurements of v sin i by 
iDelfosse et all (flQQS 'l. In GJ1156, we measure a rota- 
tional velocity three times higher than formerly reported, 
and in G J 2069 B our new value of rotational bro adening 
IS V sin I = 6kms-i while IDelfosse et al.l ()1998D reports 

V sinz — 9kms~^. We attribute the differences to higher 
data quality in our sample and to our more sophisticated 
fitting procedure. The uncertainty in v sini also depends 
on the magnetic flux with a higher uncertainty at higher 
magnetic flux. The crosstalk between magnetic flux and 
rotation velocity may be weaker in other wavelength re- 
gions (but so far this was not investigated in detail). 

Measuring magnetic flux is hampered if a target is 
rapidly rotating because resolving individual lines be- 
comes difficult in the presence of rotation. In the three 
rapid rotators Gl 493.1, LHS 3376, and GJ 1156, we 
used only the spectral regions 9946-9956 A and 9972- 
9981 A, which are particu larly well suited in rapid rota- 
tors emers fc Basrill2006' ). 

In all seven stars, we detected mean magnetic fields 
of 2kG or stronger. RB07 investigated 22 M dwarfs, 
17 of them have spectral types earlier than M7. All 
8 of them with normalized Ha luminosity larger than 
logiyHa/-^boi = —4 also have mean magnetic fields on 
the order of 2kG or more. The most active stars in 
that sample exhibit significant rotation, but the sam- 
ple contains only one star rotating more rapidly than 

V sinz = lOkms"^. Our results are in good agreement 
with the relation between mean magnetic fields and nor- 
malized Ha activity found in RB07. The (projected) ro- 
tation velocities of the three most rapidly rotating stars 
in our sample are at least a factor of two higher, but we 
see no sign of normalized Ha luminosity higher than in 
RB07, and no exceptionally high value of Bf. The main 
result of this work is that none of the three rapid rotators 
with projected rotation velocities close to 20kms~^ shows 
a mean magnetic field above 3kG. 

Can we really measure magnetic fiel ds in stars rotat- 
ing as rapidly as v sinz = 20kms~^? In IReiners fc BasTil 
|2003) we have shown that the magnetic sensitivity of 
FeH in principle allows the measurement of magnetic flux 
in stars rotating as rapidly as v sinz = 30kms~^. The 
limiting factor in rapid rotators is the achievable signal- 
to-noise ratio. In our case of "only" v sinz « 20kms~^ 
this is not a crucial problem. The differences between 
stars with strong and weak magnetic fields in the pres- 
ence of rapid rotation can be seen in the spectra we show 
in Figs. [3] and m The results of our x^ fits plotted in the 
lower panels of Fig. [3] and |4] show that in all three stars 
X^ becomes significantly larger at very low field strengths 
or if one allows for mean fields as high as 4kG. 

Our active template star Gl 873, which we are us- 
ing as comparison for our targets, has a mean magnetic 
field of about 4kG. From comparison to the spectrum 
of Gl 873 we cannot measure magnetic flux in excess of 
that value. In RB07 we show a spectrum of YZ CMi 
(M4.5, V sin z = 5 km s""'^) from which we inferred a mean 
field stronger than 4kG. This spectrum shows that in the 
presence of stronger fields the magnetically sensitive lines 
can become even more washed out following the princi- 
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pies of Zeeman broadening.* The spectrum of Gl 412 B 
(WX UMa) shows the same behavior at a rotation ve- 
locity of V sin I = 5kms~^. Thus, although we lack a 
spectrum independently calibrated to stronger magnetic 
flux to compare with, we see no way that the spectra of 
the three most rapidly rotating stars can be consistent 
with magnetic flux stronger than Bf^ 4kG. We discuss 
these stars further below. 

4.1. Saturation of the magnetic flux Bf 

It is well accepted that coronal and chromospheric ac- 
tivity saturate at high rotation rates. iMohantv fc Basrl 
((2003) showed that this relation is still valid in stars as 
late as spectral type M8.5. In mid-M stars, all stars with 
detected rotational broadening show Ha emission at the 
saturation level. In the left panel of Fig.O we plot v smi 
vs. \ogLua/Lhoi for all stars of spectral types earher 
than M7 from RB07 together with the stars from our 
new sample. This shows that the rotation-activity rela- 
tion is still valid in the combined sample. In the right 
panel of Fig.[5l we plot measured magnetic flux Bf as 
a function of v sini. If magnetic flux did not saturate, 
we would expect the values of Bf to continue growing 
with higher rotation rate. This is not observed: instead, 
the magnetic flux shows the same saturation effect as 
Ha emission. From this result we conclude that mag- 
netic flux generation does not grow any further in stars 
with saturated Ha emission. In other words, magnetic 
flux saturates in roughly the same fashion as activity, 
implying that activity does not saturate solely because 
the whole area of the star is covered with fields (/ = 1). 

4.2. Comparison to hotter stars 

How does the saturation of magnetic flux in our M- 
dwarf sample fit into the picture of rotation and activity 
in hotter stars? In order to compare activity measure- 
ments in very different stars, the Rossby number, i.e., 
the ratio of rotation period P and convective overturn 
time Tconvj is the parameter of choice. We note, however, 
that the main effect of the Rossby number is to compare 
rotation periods, and that the influe nce of the convec- 
tive o verturn time is still debatable. iKiraga fc Stepieril 
(|2007f ) provide a measurement of the rotation period 
of G1729. F or Gl 873, A D Leo, and YZ CMi, P can 
be found in ISaail (|2001h . To obtain the (projected) 
Rossby number, Ro = P/tcohv, for the other stars, we 
calculate the (projected) rotation periods for our sam- 
ple and the stars from RB07 from the projected sur- 
face rotation velocity v sinz for which we require the 
radii. Note that this only provides Ro/ smi, which is 
an upper limit of Ro. To determine the radius , we em- 
loyed the mass-luminosity relation from Delfos ^et al.l 



2000l) and the mas s -rad ius relation at an age of 5 Gyrs 
from iBaraffe et all (|1998i) . To compute th e masses, 



we used J-i nagnitudes fro n i ICutri et al.l ([2003) and dis 
tances from lHawlev et all (Il996l). Parallaxe s for G l 299 
and GJ 1286 are from [ Harrington fc DahnI ()1980D and 
lOppenheimer et"al] (|200ll ). respectively For the convec- 
tive overturn time, we adopt a value of Tca-n v = 70 d con - 
sistent with the values given for M dwarfs in lSaa^ (|200ll ). 

* We show this spectrum in Fig. 3 of RB07. In that plot, the 
ratio of the magnetically sensitive to the insensitive lines is smaller 
in YZ CMi than in our magnetically active template. 



which are taken from lGillilandl (|1986D . iKiraga fc St^piehl 
()2007D have recently calculated empirical turnover times 
for a sample of M stars. Their relation yields values 
about Tc ~ 50-100 d for the mass range considered in 
our sample. This is roughly consistent with Tc — 70 d; 
if we were using the c onvective overturn times from 
IKiraga fc St^piehl (|2007( ). the Rossby numbers would be 
about 0.15 dex larger (smaller) for the hottest (coolest) 
stars. This difference does not affect our conclusions. 
Radii, (projected) periods and Rossby numbers are given 
in TableU 

We comb ine our mea sure ments with da ta on additional 
Saad (I1996D and|Saar| (|2001h . For the stars 
19961 ) we adopt a value of Tconv = 20 d (spec- 



stars from 
from lSaarl 

tral type G0-K5) noting that the convective overturn 
time changes over this range of stars but that the differ- 
ence is not significant on the level considered here. 

The behavior of Bf with Rossby number is shown 
in Fig. [6] where we plot all the mentioned samples to- 
gether. Data from this work are plotted as filled cir- 
cles, data front RB07 as o pen c i rcles (M dwarfs), and 
data fromlHiiB (fT996l) andfSaai (l200l as crosses (note 
that crosses are hotter stars). Stars in which only up- 
per limits of V sin i were measured are not shown (non- 
detections of rotation). Fig. [6] clearly shows the depen- 
dence of magnetic flux generation with Rossby number: 
At large Rossby numbers {Ro > 0.1), i.e. in the non- 
saturated regime, magnetic flux increases with decreas- 
ing Ro (more rapid rotation). At smaller Rossby number, 
i.e., in the regime where Ha and X-ray emission saturate, 
Bf saturates as well. 

In Fig.[6l the saturated regime mainly consists of M 
dwarfs while the rising p art of the relat ion (crosses) is 
populated by hotter stars. iReinersI ()2007t ) showed that a 
drop in activity with higher Rossby number also appears 
in slowly rotating M dwarfs (his Fig. 10). A substan- 
tial piece of evidence for an intact rotation-activity rela- 
tio n among M dwarfs can be found in the data presented 
bv IKiraga fc St^piehl (|2007t ). Fr om their data, we show 
normalized X-ray activity vs. Rossby number (again as- 
suming Tconv = 70 d) in Fig.jTl Clearly, a rise of activity 
at high Rossby number {Ro — 0.1 ... 1) and a satura- 
tion plateau at lower Rossby number appears among M 
dwarfs, too. 

We conclude that the dynamo process saturates at 
Rossby number of about Ro = 0.1. No systematic in- 
crease of Bf occurs if the Rossby number is smaller by 
another order of magnitude. 

4.3. The regime of saturation 

The level of magnetic flux in the saturated regime is 
between 2 and 4kG. In two stars, we observe magnetic 
flux that may be higher than 4 kG. At least one of the two 
stars is not among those with the smallest Rossby num- 
bers in our sample. We can speculate whether the Rossby 
numbers of both stars with the highest magnetic flux are 
smaller than the Rossby numbers of the other stars in our 
sample, i.e., whether the Rossby numbers of YZ CMi and 
WX UMa are also on the order of log i?o « — 2. In the 
case of WX Uma, we only have Ro/ sinz, and for a real 
value of log i?o = — 2 the star would be observed under 
an inclination angle of z < 30°. For WX UMa, this is a 
viable option. On the other hand, a rotation period of 
P ~ 2.8 d is reported for YZ CMi, which is in good agree- 
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TABLE 3 

Masses, radii, pr,o,iected rotational velocities, periods, and Rossby numbers, magnetic flux, and normalized X-ray 
luminosities for our sample stars and the stars from rb07 (see text). 



Name 


Other 


SpT 


M/Mq 


R/Rq 


V sini 


P/ sin i 


log {Ro/ sin i) 


Bf [G] 




GJ 3379 




M3.5 


0.24 


0.25 


< 3 


>4.2 


> -1.2 


2000 


-2.86 


GJ 2069 B 


CV Cnc 


M4.0 


0.24 


0.25 


6 


2.3 


-1.5 


2500 


-2.77 


Gl 493.1 


FN Vir 


M4.5 


0.17 


0.20 


18 


0.5 


-2.1 


2100 


-3.31 


LHS 3376 




M4.5 


0.14 


0.17 


16 


0.5 


-2.1 


2000 


-3.63 


GJ 1154 A 




M5.0 


0.18 


0.20 


6 


1.7 


-1.6 


2000 


-3.28 


GJ 1156 


GL Vir 


M5.0 


0.14 


0.16 


17 


0.5 


-2.1 


2100 


-3.39 


Gl 412 B 


WX UMa 


M6.0 


0.11 


0.13 


5 


1.4 


-1.7 


> 3900 


-3.28 



Stars from RB07 



Gl 70 




M2.0 


0.35 


0.33 


< 3 


>5.6 


> 


-1.1 





< 


-4.44 


Gl 729 


V1216 Sgr 


M3.5 


0.18 


0.20 


4 


2.9*' 




-1.4 


2200 




-3.50 


Gl 873 


EV Lac 


M3.5 


0.33 


0.31 


3 


4.4*' 




-1.2 


3900 




-3.07 


Gl 388 


AD Leo 


M3.5 


0.42 


0.39 


3 


2.6'= 




-1.4 


2900 




-3.02 


Gl 876 




M4.0 


0.32 


0.31 


< 3 


>5.2 


> 


-1.1 







-5.23 


GJ 1005A 
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ment with an inclination angle close to z = 90° given the 
estimated radius and the measured rotation velocity. To 
push the value of log i?o « —1.4 to —2, either the rota- 
tion velocity must be a factor of 4 higher (i < 15° imply- 
ing that the rotation period is wrong), or the convective 
overturn time must be longer by the same a mount . Both 
options seem rather unlikely. We note that ISaaH ()200lD 
reported a magnetic flux of Bf — 3.3 kG for YZ CMi, 
a value that is somewhat lower than our result. This 
may indicate that our value does not reflect an unusu- 
ally strong average field strength in YZ CMi, but that 
the magnetic flux shows rather large scatter (either due 
to uncertainties in the measurements or temporal fluctu- 
ations). 

The easiest explanation for the two very high values 
of Bf > 4kG is that the scatter in the saturated mag- 
netic flux level is fairly large, and Bf between 2 and 5 kG 
might just be the allowed range at small Rossby numbers 
(including observational effects). We have searched for 
other parameters, such as exceptionally low gravity, that 
could cause the high flux values in the two stars. We 
did not find any particular stellar parameter that dis- 
tingiushes YZ CMi and WX Uma from other flare stars. 
In particular, age is probably of little direct importance 
for the generation of very high magnetic flux: WX Uma 
is an old disk flare star while YZ CM i is a member of the 
young disk population (|VeedeHll974f ). 

Even if the two strongest magnetic flux measurements 
are due to exceptionally small Rossby numbers, this 
would not explain the saturation of Bf between log Ro = 
—2 and —1. Across an order of magnitude in Ro, Bf 



varies by at most a factor of a few. This is in strik- 
ing contrast to the nearly hundredfold increase in Bf in 
going from Ro « 1 to Ro « 0.1. We have concluded 
that both magnetic flux and chromospheric Ha (as well 
as coronal X-ray) emission saturate at small Ro. With 
the large uncertainties in Bf, we cannot reliably deter- 
mine whether there is any super-saturation effect, with 
magnetic flux declining in the most rapid rotators. 

Although both emission and magnetic flux appear to 
saturate at small Ro, we observe a large scatter of mag- 
netic flux among the stars with small Ro, and likewise 
some scatter in Ha emission. We thus examined whether 
a relation between Ha emission and magnetic flux still 
exists at very small Rossby numbers. In Fig.[8l we plot 
normalized Ha emission vs. magnetic flux for our stars. 
As expected, logLnaZ-^^boi grows as Bf grows from zero 
to 2 kG. Beyond that, i.e. in the regime of saturated mag- 
netic flux, no further increase in chromospheric emission 
is observed, although Bf grows as large as 4kG. No ob- 
vious correlation exists between chromospheric emission 
and magnetic flux in the regime of saturation. This could 
mean that saturation of chromospheric emission occurs 
above a certain level of integrated flux on the star (for 
example, because the filling factor reached unity). The 
field could then grow stronger without causing further 
heating. This is in fact not ruled out by our poorly con- 
strained understanding of chromospheric heating mech- 
anisms. 

As of now, the available data do not permit us to dis- 
tinguish reliably between these two possibilities. But, as 
noted above, the most rapid rotators in our sample prob- 
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ably do not possess magnetic fluxes higher than among 
other stars in the saturated regime: in Fig.[8l stars with 
the smallest Rossby numbers lie at the \aw-Bf end of 
the plateau. Neither are they the stars with the high- 
est normalized chromspheric emission. Thus the lack of 
tight correlation between emission and Bf in the satu- 
ration regime should not be construed as evidence that 
magnetic flux continues to be correlated with rotation, 
but decouples from emission. Instead, our data suggest 
that any variation of Ha emission or Bf in the saturation 
regime must either be random - i.e. due to the obser- 
vational scatter - or depend on parameters other than 
rotation. At this stage, we note only that the plateau in 
Fig. [H] is consistent with scatter alone. 

5. SUMMARY 

Using absorption lines of molecular FeH we have mea- 
sured the magnetic flux in seven M stars that are known 
X-ray sources. With one exception they were known as 
rapid rotators as well. We reinvestigated the projected 
rotational velocities v sini and measured chromospheric 
emission in Ha. All stars proved to be strong Ha emit- 
ters. In our analysis of v sini, we found some inconsis- 
tencies with former literature, which we ascribe to our 
more sophisticated (direct) fitting method and better 
data quality. 

All of our target stars show strong Ha emission at the 
saturation level (logLna/iboi > —4). We detected mag- 
netic fields of kG-strength in all our targets. While in 
less rapidly rotating stars. Ha and magnetic flux corre- 
late with rotation rate, no such correlation is observed in 
our sample. 

In contrast to sun-like stars of spectral types G and 
K, M dwarfs rotating at v sini « lOkms"^ have very 
small Rossby numbers. Thus, our sample targets add 
to the amount of stars with measured magnetic flux in 
the regime of small Rossby number, i.e. in the regime 
of saturated magnetic activity. Our primary goal was 
to determine whether magnetic flux in this regime con- 
tinues to grow with more rapid rotation, as sometimes 
suggested (e.g.. ISaail 119961 ) . or instead saturates in the 
same manner as coronal and chromospheric emission. 

Our main conclusion is that around a Rossby number 
of Ro ~ 0.1, magnetic flux saturates at approximately 
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Bf — 3kG. Below Ro ^ 0.1, Bf does not grow stronger 
with decreasing Rossby number. In looking at the ef- 
fect of Rossby number {Ro = P/rconv), we are primarily 
sensitive to the effects of rotation, because the sampled 
rotations vary by a factor of more than ten, while the 
convective overturn time probably changes less than a 
factor of two in the range of M stars we are considering. 

In the regime of saturated magnetic flux and chromo- 
spheric emission, we still observe a strong scatter in mag- 
netic flux. The interpretation of this feature is not clear. 
It may be that below the critical Rossby number, normal- 
ized Ha emission is not sensitive to changes in Bf . Alter- 
natively, in light of the large systematic errors of our Bf 
measurements, the scatter of Bf in the saturated regime 
may be fully explained by observational uncertainties. 

Our results indicate that the strengths of stellar mag- 
netic fields, and not merely their filling fractions on the 
surface alone, reach a maximum at a certain rotation 
rate. The scatter in magnetic flux among the stars 
with very small Rossby numbers, however, is substantial. 
Typical values are between 2 and 4kG. Two stars show 
magnetic flux stronger than 4kG, and we cannot exclude 
that individual stars can generate magnetic flux stronger 
than that. Determining what sets the maximum value - 
whether equipartition with either the atmospheric pres- 
sure or the turbulent velocity field, or some more subtle 
effect - remains a challenge for future observations and 
theories. 
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Fig. 1. — Top panel: Data and fit of the two slow rotators GJ 3379 (left) and GJ 2069B (right). The data are shown in black. Our fit 
to the data for the case of no magnetic flux is overplotted in blue, very strong magnetic flux (B/ = 3.9 fcG) in red, and the best fit with 



intermediate flux values in green. Bottom panel: x l^'iidscapes showing the goodness of flt as a function of v sini and Bf . 
color indicates good fit quality, bright/yellow color means bad fit. The white contour marks the formal 3o" region, i.e., < 
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Fig. 2. — Date and fit (top panel) and x^-landscapes (bottom panel) as in Fig.[T]for the two intermediate rotators GJ 1154A (left) and 
GJ 412B (right). 
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Fig. 4.— As Fig.|3]for the rapid rotator LHS 3376. 
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Fig. 5. — Left panel: Normalized Ha activity as a function of v sin j in M stars; right panel: Magnetic flux Bf as a function of v sini. 
Filled circles are from this work, triangles come from RB07. Downward arrows indicate upper limits, numbers give the number of multiple 
measurements with the same results. The two lower limits of Bf are indicated with upward arrows. 
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Fig. 6. — Magnetic flux Bf as a function of Rossby number. Data are from lSaail lil9M ). [Saa^ ll200ll ). RB07, and from this work. 11 
stars without measured rotation periods and with no detection of rotational broadening (v sini < 3kms~^), i.e. lower limits of Ro, are 
not shown. They would form a vertical line at about log i?o = — 1 but probably lie on top of the rising part of the c orrela t ion (s ee IReiner? 

t7f). Data from this work are plotted as filled circles, data from RB07 as open circles (M stars), and data from [Saail ll 19961 ) and ISaa" 
31) as crosses (spectral types G0-M2). Rotation rates from this work and from RB07 are calculated from v sini implying that open 
and filled circles are Ro/ sini hence upper limits of Ro. 
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Fig. 7. — Normalized X-ray activity as a function of Rossby number in M stars (using Tconv 

= 70 d) from lKiraea fc StepienI i200% . This 
plot is essentially the same as their Fig. 7 but in logarithmic units so that the saturation plateau becomes clear. 
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Fig. 8. — The correlation between Ho activity and magnetic flux Bf in M stars. Symbol size scales with inverse Rossby number (large 
symbols have small Ro). Triangles are from RB07, circles are from this work. 



